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Approaches  to  Computer  Modeling  of  Phosphate  Hide-Out 


X.  A.  S.  Hardy  and  J.  C.  Cooper 
Cheaistry  Division 
Naval  Research  Laboratory 
Washington,  DC  20375 


INTRODUCTION 

Phosphate  is  nsed  as  a  vater  treatment  chemical  in  navy 
boilers.  The  main  purpose  of  this  treatment  is  to  remove  calcium 
from  solution  as  soft  nonadherent  calcium  phosphate  sludge  rather 
than  as  hard  scale  forming  species  such  as  calcium  carbonate, 
calcium  sulfate,  or  calcium  silioate.  The  calcium  phosphate 
formed  can  then  be  removed  from  the  boiler  by  blowdown.  Addition 
of  NalP04  and  NasHP04,  called  the  coordinated  phosphate  or  cophos 
treatment,  also  affords  pH  control  as  the  phosphate  acts  as  a 
buffer  to  keep  pH  at  a  value  above  which  acid  corrosion  occurs 
and  below  which  caustic  corrosion  becomes  significant. 

Difficulties  are  sometimes  encountered  in  the  use  of  this 
cophos  treatment,  and  phosphate  levels  decrease  below  the  desired 
values.  Some  of  these  decreases  are  due  to  contamination  by 
shore  or  sea  water,  by  outward  leakage,  or  by  carryover. 
Sometimes,  however,  decreased  phosphate  levels  cannot  be 
explained  by  these  causes;  this  is  called  phosphate  or  chemical 
hide-out  and  it  is  characterized  by  a  decrease  in  phosphate 
concentration  at  higher  firing  rates.  The  hide-out  can  be 
reversed  by  lowering  the  firing  rate.  This  makes  control  of  the 
phosphate  levels  very  difficult  because  addition  of  phosphate  at 
the  higher  firing  rates  to  keep  the  phosphate  levels  within  the 
established  limits,  causes  the  phosphste  concentrations  to 
increase  beyond  these  limits  upon  decreasing  the  firing  rate.  In 
addition,  pH  control  is  lost  and  pH  rises  into  the  caustic 
region. 

(taputer  modeling  of  phosphate  hide-out  in  boiler  water 
could  predict  when  the  phenomenon  will  occur  and  suggest 
Manuscript  approved  March  28, 1984. 
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corrective  actions  for  prevention.  The  accuracy  and  predictive 
value  of  the  eoaputer  model  depends  on  the  quality  of  the  data 
and  the  validity  of  the  aatheaatical  alogrithas  used  in  the 
eoaputer  prograa.  Developaent  of  such  a  aodel  requires,  first, 
deteraination  of  the  significant  cheaieal  species,  inforaation 
about  their  oonoentrations,  and  an  assessaent  of  available  high 
teaperature  kinetic  and  theraodynaaie  data  for  these  species. 

Given  sufficient  data,  a  aodel ing  aethod  can  then  be  selected  and 
complex  software  developed. 

Suaaarized  in  this  report  is  such  a  survey  of  available 
inforaation.  After  reviewing  the  current  understanding  of 
phosphate  hide-out,  data  on  species  and  species  concentrations 
were  collected.  High  teaperature  theraodynaaie  data  were 
compiled.  Available  eoaputer  programs,  useful  for  aodeling  aulti- 
species,  high  teaperature  solutions,  were  obtained  and  modified. 
The  sufficiency  of  the  data,  and  the  adequacy  of  the  programs 
were  evaluated.  Initial  attempts  made  to  develop  a  basic  aodel 
using  these  data  and  programs  are  presented.  Keeoaaendations  are 
aade  concerning  collection  of  iaportant  data  and  further 
developaent  of  programs.  The  review  of  the  current  understanding 
of  phosphate  hide-out  follows. 


BACKGROUND  ON  PHOSPHATE  HIDE-OUT 

Several  explanationa  are  reported  for  phosphate  hide-out. 

It  has  been  suggested  [1-5]  that  a  fora  (or  forms)  of  sodium 
phosphate  exists  which  is  less  soluble  at  higher  temperatures  and 
that  this  causes  phosphate  to  hide  out  at  high  temperatures  and 
to  redissolve  at  lower  temperatures.  Interaction  of  iron  and 
phoaphate  as  NaFeP04  [4-7]  and/or  Fes(P04),  [1,7]  has  also  been 
suggested  to  explain  the  loss  of  phosphate  froa  solution.  Hide¬ 
out  has  been  seen  to  occur  to  a  greater  degree  in  recently 
cleaned  boilers  which  were  not  sufficiently  passivated  [4].  Hide¬ 
out  alao  occurs  when  waterside  deposits  are  present  [8] , 
especially  when  these  inolude  such  species  as  Fe,  Cu,  Ni,  and  Zn 
[4]  and  in  areas  where  corrosion  and  concentrating  filas  are 
found  [9]  such  as  crevices,  pits  and  bends  [8] .  A  better 
explanation  of  the  hide-out  aay  be  obtained  froa  ooaputer 
aodeling  with  the  proper  paraaeters.  Such  a  aodel  could  be  used 
to  control  and  predict  this  unwanted  phenoaenon. 

Several  reports  describe  attempts  to  develop  aodel s  of  steam 
generating  systeas  [10-13]  and  of  industrial  cooling  water 
systeas  [14-17] .  Purposes  and  orientations  vary  but  the  chief 
goal  of  all  such  aodel s  is  to  reduce  the  relationships  of  the 
species  in  a  solution  to  aatheaatical  terms,  or  equations,  which 
can  be  used  to  siaulate  the  operation  and  interrelationships  of 
these  species.  These  siaulations  esn  then  be  used  as  an  approaoh 
to  understand  such  things  as  corrosion  [12]  and  scale  foraation 
[14] ,  and  to  predict  the  changes  which  will  be  caused  by  such 
perturbations  as  the  addition  of  cheaieal s  and  changes  in 
teaperature.  This  inforaation  can  then  be  used  to  choose  safe 


operating  limits  and  appropriate  water  treatment.  None  of  these 
reported  attempts  has  focused  on  phosphate  hide-oat. 


SPECIAXION,  CONCENTRATIONS,  AND  CONDITIONS 

The  first  step  necessary  in  preparing  a  aodel  is  to  iteaize 
the  species  in  solntion,  especially  those  which  will  or  aay  be 
iaportant  in  phosphate  hide-ont.  Sources  of  this  inforaation  are 
disevssed  in  Appendix  A  and  a  list  of  these  species  is  available 
in  Table  A-l  of  that  Appendix.  These  are  individual  species  and 
theraocheaieal  data  are  needed  for  the  ooabinations  and  reactions 
of  thea. 

Besides  inforaation  on  the  types  of  species  found  in  a 
boiler,  inforaation  is  also  needed  on  the  concentrations,  or 
ranges  of  concentrations,  of  those  species.  On  the  whole  there 
is  not  aueh  inforaation  available  to  us  systeaatically 
suaaarizing  the  types  and  aaounts  of  various  species  in  boiler 
water.  Appendix  A  suaaarizes  the  sources  of  such  inforaation  and 
the  data  which  are  available.  The  quantitative  data  froa  these 
sources  are  presented  in  Table  A- 2  of  that  Appendix. 

Although  this  inforaation  provides  a  starting  place,  it  aust 
be  stated  that  a  aore  comprehensive  analysis  of  the  species  found 
in  boiler  water  while  operations  are  noraal  and  while  phosphate 
hide-out  is  occuring  is  desirable. 

The  conditions  found  in  the  boilers  included  a  range  of 
temperatures  from  25*C  to  300*C.  If  conductivity  is  kept  within 
the  suggested  levels  (less  than  300  pahos/ca) ,  then  ionic 
strength  is  low  (probably  less  than  0.01  M) .  However,  fault 
conditions  such  as  inleakage  or  concentrating  filas  can  raise 
ionic  strength  significantly. 

KINETIC  AND  THERVOCHEHICAL  DATA 

For  a  coaplete  aodel  of  any  solution  both  kinetic  and 
equilibriua  thermodynamic  data  are  needed.  It  is  important  to 
include  both  approaches.  For  example,  one  reaction  aay  be 
energetically  favored  but  will  still  not  occur  within  the  given 
tiae  frsae  because  the  process  is  very  slow.  Such  a  case  is  the 
precipitation  of  Ca,(P04)l0H,  hydroxyapatite,  which  is  the 
thermodynamically  aost  stable  calciua  phosphate.  However, 
depending  on  the  pH,  it  is  not  always  the  first  form  to 
precipitate  because  kinetically  it  is  a  slower  reaction. 

Nanoollas  [18]  has  studied  the  kinetics  of  these  precipitations 
and  has  discussed  the  interplay  of  kinetics  snd  thermodynamics 
for  several  calciua  phosphate  foras.  However,  the  work  was  done 
at  low  temperature  (37*0  and  data  for  high  temperatures  are  not 
available.  Other  data,  also  available  only  at  low  temperatures, 
include  that  on  rust  formation  [19] .  The  rate  data  for  coaplex 
and  ion-pair  formation  and  on  the  dissolution  and  precipitation 
reactions  at  high  temperatures  are  not  available. 


This  doss  not  assn  that  Modeling  is  iapossible.  Other 
aodsls  based  on  theraodynaaics  have  been  successful  in  asking 
useful  predictions  for  boiler  water,  oooling  water,  and  natural 
waters  [20].  It  is  also  likely  that  at  high  teaperatures  aost 
kinetics  are  sufficiently  fast  that  good  aodels  can  be  developed 
using  a  theraodynaaic  approach.  Te  decided  to  continue  to 
evaluate  the  availability  of  the  necessary  high  teaperature 
theraodynaaic  data  and  to  pursue  Modeling  froa  this  approach. 

A  large  body  of  theraocheaical  data  at  298*K  (25*0  has  been 
coapiled.  Data  for  the  high  teaperatures  found  in  boilers  (up  to 
300*0,  however,  are  spread  throughout  the  literature.  Vater 
itself  behaves  very  differently  at  high  teaperatures.  It  is, 
therefore,  necessary  to  disouss  the  nature  of  vater  and  of 
solutions  at  high  teaperature,  and  to  suaaarize  aethods  of 
dealing  with  various  foras  of  available  data  before  describing 
the  collected  data.  These  are  discussed  in  Appendix  B. 

To  aodel  chemistry  in  boilers  the  availability  of  high 
teaperature  data  for  the  species  of  interest  in  boilers  aust  be 
assessed.  A  ooapilation  of  available  data  is  presented  in 
Appendix  C;  it  includes  sources  of  this  theraocheaeial  data,  a 
description  of  the  foras  in  which  these  sources  present  data,  and 
Tables  C-l  through  C-9  in  which  the  actual  data  are  organised  by 
species  type. 

COMPUTER  PROGRAMS  AND  MODELING 

Having  snuaarixed  the  available  data,  a  consideration  of 
available  prograas  and  other  aodel ing  approaches  led  to  the 
selection  of  two  prograas  with  which  to  begin  initial  Modeling. 
The  program  COMICS,  developed  by  Perrin  and  Sayoe  [20]  and 
revised  by  Cooper  [21]  was  already  available  in  this  laboratory. 
The  program  calculates  the  concentrations  of  free  metals,  free 
ligands  and  of  the  coaplex  species  in  solution,  including  aixed 
species,  such  as  MM'L,  hydrolyzed  species,  such  as  M(OH)s, 
protonated  species  such  as  NHL,  and  polynuclear  species  such  as 
MaL, ,  where  M  and  M'  represent  different  aetals  and  L  represents 
soae  ligand.  Using  foraation  constants,  total  ligand  and  aetal 
concentrations,  end  pH,  the  prograa  calculates  the  concentrations 
of  up  to  forty  species  derived  froa  four  different  aetal  ions  and 
four  ligands,  and  will  plot  these  results.  The  prograa  aakes  no 
corrections  for  high  teaperature  except  of  course  that  the 
foraation  constants  used  oan  be  the  high  teaperature  values  since 
they  are  input  by  the  user.  Also,  no  precipitation  reactions  can 
be  included.  The  second  program  chosen  for  use  is  WATEQF  [22] 
which  is  a  Fortran  prograa  written  for  the  calculation  of  the 
cheaical  equilibria  in  natural  waters;  it  was  developed  for  use 
by  the  United  States  Geological  Survey  (USGS) .  Included  in  the 
prograa  are  193  reactions,  both  aqueous  foraation  and 
precipitation  reactions.  The  prograa  includes  corrections  for 
higher  teaperatures  in  the  calculation  of  activities  froa  ionic 
strength,  soae  high  teaperature  equilibia  data,  and  allowance  for 
the  change  in  concentrations  due  to  the  chsnge  in  density  as 


temperature  changes. 

Becanse  COMICS  was  available  and  had  the  capacity  to  handle 
so  many  species  in  solution,  initial  calculations  were  done  with 
it.  The  first  calculations  were  aade  of  the  ionization  of 
phosphoric  acid  at  low  (25*0  and  high  (300*C)  tenperature. 
Foraation  oonstant  data  were  first  calculated  using  the  entropy 
and  enthalpy  data  [23]  in  Table  C-3  using  equations  (B-3  and  B-4) 
and  assuaing  a  zero  heat  capacity.  It  is  instructive  to  coapare 
the  log  of  the  equilibrium  constants  so  calculated  with  those 
aeasured  by  Mesmer  and  Baes  [24,25]  for  reactions  (1-3): 


H+ 

+ 

h,po4- 

-  h,po4 

(1) 

fit 

+ 

hfo4»- 

-  h,fo4- 

(2) 

Bf 

+ 

P04»- 

-  hpo4»- 

(3) 

For  reaction  (1)  the  calculated  value  at  25*C  is  2.12,  while  the 
aeasured  value  is  2.145.  At  300*C  the  calculated  value  is  2.78 
while  the  aeasured  constant  is  4.40.  For  reaction  (2)  at  25*C 
the  calculated  value  is  7.21  and  the  aeasured  value  is  7.197 
while  at  300°C  the  calculated  oonstant  is  5.93  while  the  aeasured 
is  8.945.  Note  that  aeasured  and  calculated  values  are  very 
close  at  25*C  but  are  not  so  at  300*C.  The  difference  clearly 
deaonstrates  that  heat  capacity  aust  be  included  at  the  higher 
temperatures,  and  that  there  is  a  greater  degree  of  ion 
association  at  the  higher  temperatures.  The  value  of  the  log  of 
the  third  constant,  reaction  (3),  is  not  given  at  300°C  because 
it  was  too  difficult  to  extrapolate  to  solutions  of  low  ionic 
strength.  For  the  calculation  with  the  COMICS  program  a  value  of 
12.7  was  assumed  for  reaction  at  300*C  as  this  represents  a  good 
approximation  based  on  the  data  for  150*C  [24] .  The  total 
concentration  of  phosphate  was  taken  as  25  ppm  or  2.5  x  10-*  M. 
Figures  1  and  2  summarize  the  calculations  done  at  25*C  and 
300*C.  A  comparison  of  the  two  plots  indicates  that  at  the 
higher  temperature  the  phosphate  stays  bound  to  the  proton  up  to 
higher  pH  values.  This  means  that  in  the  pH  range  of  normally 
operating  boilers  the  phosphate  is  found  as  HP04*-  or  HaP04-. 

This  is  consistent  with  the  decreased  solvating  capacity  of  water 
at  high  temperatures,  such  that  a  higher  degree  of  ion 
association  and  complex  formation  is  seen. 

Conclusions  which  can  be  drawn  from  this  COMICS  calculation 
with  the  ion  association  constants  for  phosphate  are:  heat 
capacity  must  be  included  in  the  calculation  of  high  temperature 
formation  constants  and  at  higher  temperatures  there  are  higher 
concentrations  of  the  protonated  forms  of  phosphate. 

COMICS  was  also  used  for  the  calculation  of  metal  bound 
phosphate  using  the  constants  reported  for  25*C  by  Sillen  and 
Nartell  [25,23]  in  Table  C-3.  The  concentrations  of  dissolved 
metal  ions  were  chosen  based  on  reported  vslues  [20]  shown  in 
Table  A-2.  These  values  are:  Ca,  .03  ppm  -  .7  x  10-*  M; 


Fe,  .09  ppm  -  1.6  x  10-*  M;  Co,  .17  ppm  -  2 .7  x  10-*  M;  Mg, 

.18  ppa  ■  7.4  x  10-‘  M.  A  separate  calculation  was  done  with 
each  aetal  and  25  ppa  phosphate.  Another  calculation  was  then 
done  for  the  aixtnre  of  all  the  aetals  with  the  phosphate  nsing 
the  saae  concentrations  described  above.  The  acid  association 
values  used  were  those  of  Martell  and  Sillen  [23] .  Figure  3 
plots  the  calculated  percentages  of  species  for  the  calculation 
done  with  ealeiua  and  phosphate.  The  species  HaP04-,  HP04*-,  and 
P04 *-  were  calculated  as  per  cent  total  phosphate  whereas  free 
Ca*+  and  CaHF04  were  calculated  and  plotted  as  per  cent  total 
calciua.  Because  phosphate  is  in  suoh  large  excess  (26  fold) 
over  the  ealeiua,  the  aaonnt  of  calciua  hydrogen  phosphate 
present  is  very  saall  and  represents  less  than  1%  of  the  total 
phosphate.  The  curves  for  the  phosphate  remain  unchanged  from 
that  calculated  without  metal  present  as  shown  in  Figure  1.  The 
curves  for  the  metal  speoies  in  Figure  3  indicate  that  most  (90%) 
of  the  calciua  reaains  unbound  to  the  phosphate.  This  holds  true 
for  aagnesiua  and  copper  also.  However,  the  calculation  with 
iron  shows  that  uost  090%)  of  the  iron  is  bound  as  FeHF04. 

The  results  for  the  aixtnre  of  all  the  metals  with  phosphate 
are  plotted  in  Figures  4  and  5.  For  Figure  4  the  results  are 
calculated  as  per  cent  total  phosphate.  Again  the  curves  reaain 
essentially  unchanged  froa  Figure  1.  The  aaonnt  of  phosphate 
bound  to  the  aetals  reaains  saall  even  for  the  mixture.  This  is 
due  to  the  low  concentration  of  the  aetals  and  to  the  fact  that 
only  iron  binds  phosphate  very  strongly.  Figure  5  plots  the 
results  for  each  aetal  species  in  the  mixture  calculated  as  a  per 
cent  of  that  total  metal  concentration.  The  conclusion  drawn 
froa  this  is  the  saae  as  that  for  the  individual  calculations; 
all  of  the  aetals  are  essentially  free  except  iron.  One  further 
calculation  was  done  with  the  species  concentrations  coaputed  for 
calciua,  phosphate,  and  pH.  The  most  stable  precipitate  of 
calciua  phosphate  is  Ca f (P04 ) , (OH)  and  the  solubility  product  for 
this  at  25*C  is  (4): 

K,p  -  [Ca]*[P04]»[0H]  -  10-**  (4) 

Using  the  coaputed  concentrations  at  pH  10,  with  only  calcium  and 
phosphate  present  in  solution,  the  above  product  was  found  to  be 
1.1  x  10-**.  When  the  other  aetals  were  present  the  calculated 
product  is  1.1  x  10-**.  Thus  precipitation  will  occur  under 
these  conditions  and  will  be  exaserbated  by  the  competitive 
binding  of  phosphate  by  other  aetals.  This  emphasizes  the  need 
to  include  precipitation  reactions;  COMICS  does  not  have  this 
capability.  Also  only  low  temperature  data  are  available  for 
these  aetal  phosphate  species.  It  is  expected,  based  on  the 
nature  of  water  at  high  temperature,  that  the  degree  of  metal- 
hydrogen  phosphate  association  will  increase  with  increasing 
temperature.  If  corrections  were  also  made  for  the  higher 
concentrations  of  species  found  in  orevices  and  below  porous 
deposits  then  this  iron  hydrogen  phosphate  association  could 
prove  very  important  in  phosphate  hide-out. 


Inasmuch  at  COMICS  does  not  include  precipitation  reactions, 
present  effort  involves  calculations  with  VATEQF.  Although 
VATEQF  includes  many  features  necessary  to  the  modeling  of  a 
system  at  high  temperature  as  was  mentioned  before,  some 
adaptations  are  required.  The  program  does  not  include  all  of 
the  high  temperature  data  needed  and  these  must  be  added. 
Modifications  of  this  program  for  available  computer  systems  and 
for  applicability  to  high  temperatures  are  in  progress.  It  is 
expected  that  this  program  will  offer  a  great  deal  of  capability 
to  the  modeling  of  boiler  water. 

SUMMARY 

The  purpose  of  this  study  was  to  evaluate  the  feasibility  of 
modeling  phosphate  hide-out  in  navy  boilers.  It  was  necessary  to 
collect  the  available  appropriate  data,  to  obtain  and  modify 
computer  programs  which  can  accurately  model  multi-oomponent 
equilibria  at  high  temperatures,  and  to  evaluate  the  sufficiency 
of  the  data  and  the  adequacy  of  the  programs. 

Data  on  speciation  and  the  concentrations  of  species  are 
summarised  in  Tables  A-l  and  A-2.  These  data  are  probably 
sufficient,  though  a  systematic  and  thorough  analytical 
evaluation  would  increase  the  data  base  and  insure  the 
reliability  of  the  calculated  results.  Such  a  study  should 
evaluate  the  species  and  their  concentrations  in  boiler  water 
under  normal  operating  conditions  and  under  the  condition  of 
phosphate  hide-out. 

Very  few  applicable  high  temperature  kinetic  data  are 
available.  This  does  not  mean  that  modeling  is  impossible  as 
many  systems  have  been  successfully  modeled  using  only  a 
thermodynamic  approaoh.  The  available  high  temperature 
thermochemical  data  are  summarized  in  the  Tables  of  Appendix  C. 
The  data  which  are  not  available  and  which  msy  be  of  significance 
are  the  formation  constants  and  solubility  products  for  metal 
phosphate  species.  The  compounds  which  may  be  of  special 
importance  are  iron  phosphates,  sodium  iron  phosphate,  and 
calcium  phosphates  at  high  temperature. 

Two  computer  programs  have  been  used  that  calculate 
concentrations  of  species  in  multi- component  equilibria.  Initial 
models  were  done  with  COMICS  which  does  not  include  solubility 
products.  Also,  high  temperature  corrections  such  ss  one  for 
ionic  strength  are  not  included  in  the  program.  VATEQF  is  still 
being  modified  for  this  problem,  for  high  temperatures,  and  for 
the  HP  1000  and  VAX  11/750  computer  systems.  VATEQF  has  all  of 
the  necessary  components  for  modeling  boiler  water,  or  can  be 
modified. 

Initial  modeling  done  with  COMICS  led  to  several 
conclusions.  First,  calculation  of  high  temperature  formation 
constants  using  heat  capacity  information  is  essential.  Second, 
there  is  a  greater  tendency  for  phosphate  to  be  protonated  at 


higher  temperatures.  Proa  the  modeling  done  vith  metal  ions  and 
phosphate  at  25*C,  it  can  he  concluded  that  iron  has  a  much 
higher  affinity  for  hydrogen  phosphate  than  does  calcium, 
magnesium,  or  copper.  This  affinity  is  expected  to  be  even 
higher  at  high  temperatures,  although  such  thermocheaical  data 
are  not  available.  The  inclusion  of  solubilty  data  in  these 
calculations,  not  presently  possible  vith  the  program  COMICS,  is 
essential  in  a  calculation  vith  the  large  number  of  species 
expected  for  an  accurate  and  predictive  model  of  boiler  vater 
chemistry.  Based  on  a  rough  calculation,  it  is  expected  that  the 
presence  of  other  metals  vill  increase  the  tendency  of  calcium 
phosphate  to  precipitate.  Data  and  computer  programs  have  thus 
been  collected  and  evaluated,  and  some  initial  modeling  has  been 
accomplished. 
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Fig.3-  Spociation  of  calcium  and  phosphata  a*  a  function  of  pH 

(0,H2P04*;  ♦,HP042*j#,  PO^jthaca  thraa  spacias  ora  calculatad 
as  par  cant  total  phosphata;  «,  CaHP04;  •*  ,Ca2*;  thasa  two  spacias 
ara  calculatad  as  par  cant  total  calcium  ) 


Flp.4  -  Spoc  lotion  of  calcium,  magntilum,  Iron,  coppor,  and  phoiphat 
as  a  function  of  pH 

(0,H,P04~:  +,HP042\  «,P043',*  sum  of  MHP04  whom 
M  is  Ca,Mo,Fa,  and  Cu) 


PER  CENT  SPECIES 


1 


xxxxx 

S09oo 

4»  U 


*  "’'■►°0000000,ff&cj000000000 


.'^♦;ooooooo°?iS.s.ooooooooo0;»»* 

S°  I  frH  F.’*  i  i  °6 


3 


APPENDIX  A 


SPECIATION  AND  SPECIES  CONCENTRATIONS 


Essential  to  the  development  of  a  good  nodal  is  a 
coapilation  of  the  speeies  found  in  boiler  water.  Fron  the 
enggestione  of  possible  oanses  for  phosphate  hide-out  described 
in  the  seetion  on  BACKGROUND,  the  species  phosphate,  iron, 
sodinai,  eopper,  nickel,  and  sine  are  considered  to  be  important. 

A  Naval  Ship’s  Teohnieal  Manual  (NSTM)  [8]  lists  other 
contaminant  metals  and  anions  found  in  source  water  and  in  the 
analysis  of  sludge.  Other  species  were  discussed  at  a  workshop 
of  the  American  Society  of  Mechanical  Engineers  (ASME)  [27] .  In 
other  reported  models  [12,13,28]  other  contaminants  are 
mentioned.  Two  reports  [29,30]  include  still  other  species  found 
in  navy  boilers.  A  list  of  all  of  these  species  is  given  in 
Table  A-l. 

Also  essential  is  information  on  the  concentrations  of  these 
species.  The  same  sources  with  qualitative  information 
[8,29,30,27]  as  well  as  ships'  logs  provide  some  quantitative 
data.  These  data  are  summarised  in  Table  A-2. 

Additional  qualitative  and  quantitative  data  on  boiler  scale 
is  useful.  Deposits  on  boilers  were  analyzed  and  found  to 
include  Fe,  Si,  Ca,  Mg,  S,  Pb,  Cu,  Zn,  Ni,  and  P  [29].  The 
quantitative  analysis  of  this  deposit  showed  that  0.2%  was  Pb  and 
the  authors  suggest  that  the  probable  components  of  the  deposit 
are  as  follows:  NaaS04,  99  ppm;  NaCl,  146  ppm;  I,C0,,  36  ppm; 
KaS04,  13  ppm;  SiOa,  5  ppm;  NaOH,  37  ppm;  HCO,-,  present.  The 
NSTM  also  defines  the  individual  contribution  made  to  boiler 
water  conductivity  by  the  species  C1-,  S04*-,  H00,-,  Na+,  Ca*+, 
Mg*+,  and  OH-  [8] . 


TABLE  A-l 


SPECIES  PRESENT  IN  BOILER  WATER  • 

ANIONS 

CATIONS 

NEUTRAL 

P04»- 

Bf 

o» 

co,»- 

Ca*+ 

CO, 

Cl- 

Fe*+.Fe»+ 

N, 

Si04«- 

Mg»+ 

SiO, 

SiO  *- 

Cn*+»Cn+ 

(NH,) 

0H- 

Na+ 

<B,) 

SO  *- 

Zn*+ 

0*- 

Ni»+ 

HO,- 

Mn*+ 

(Cr*+) 

(Ti»+) 

( Sn*+) 

(Pb>+) 

(Al*+) 

<K+) 

Those  in  parentheses  are  expected  to  be  less  important. 


APPENDIX  B 


DISCUSSION  OF  WATER  AND  SOLUTIONS  AT  HIGH  TEMPERATURE 
AND  SUMMARY  OF  METHODS  FOR  DEALING  WITH  VARIOUS  FORMS 
OF  HIGH  TEMPERATURE  THERMOCHEMICAL  DATA 


WATER 

The  highly  structured  nature  of  water  at  rooa  temperature 
due  to  hydrogen  bonding  is  lost  at  higher  temperatures  as  the 
kinetic  energy  increases.  This  change  in  the  degree  of  structure 
causes  dramatic  changes  in  several  properties.  Density  decreases 
(.9971  at  2 5*C,  .7125  at  300*0.  Visoosity  decreases.  The 
dielectric  constant  decreases  (78.381  at  25*C,  20.052  at  300°C) 
showing  a  change  from  a  very  polar  medium  to  one  that  is  similar 
to  non-polar  liquids.  The  solvating  nature  of  water  thus 
changes.  At  25°C  water  hydrates  ions  and  polar  molecules  easily, 
fitting  them  into  the  structure.  At  higher  temperature  water  is 
not  as  good  at  hydrating  species  and  there  is  not  as  much 
structure  for  them  to  fit  into.  Electrolytes  thus  become 
weaker.  This  affects  the  thermodynamic  functions  of  solute 
species  and  this  change  is  reflected  in  the  heat  capacity 
function.  A  similar  change  in  heat  capacity  is  not  seen  for 
homogeneous  solutions  in  non-polar  solvents.  For  the  aqueous 
solutions  found  in  boilers,  these  changes  must  be  included  and 
heat  capacity  data  at  high  temperatures  must  be  available.  It  is 
essential  that  conclusions  not  be  based  on  low  temperature  data. 

The  equilibrium  for  the  ionization  of  water  has  been  well 
studied  at  high  temperatures: 

H,0  -  H+  +  OH"  (B-l) 

Measured  pX^  values  are  available  up  to  350°C.  The  value  of  K 
decreases  as  temperature  increases,  demonstrating  the  change  or 
water  from  a  very  polar  and  ionizable  solvent  to  one  which  is 
less  so.  This  change  in  K  also  means  that  the  values  of  pH  also 
vary  greatly  with  temperature  and  modeling  must  include  this. 


IONIC  STRENGTH 

The  effect  of  ionic  strength  on  solution  equilibria 
increases  with  increasing  temperature  due  to  changes  in  the 
nature  of  water.  One  way  to  include  this  in  calculations  is  to 
use  an  extended  form  of  the  Debye-Huckel  equation  (B-2)  to 
evaluate  the  activity  coefficient,^,  which  measures  the  non¬ 
ideality  of  the  electrolyte:  ~ 

log  T+  -  ~  A|Z+Z-|*I1/2/(l  +  a  •  I1/2)  +  B  •  m  (B-2) 

Values  for  A,  the  Debye-Huckel  limiting  slope,  increase  with 


temperature  and  are  known  [25,31]  .  The  values  of  Z  are  the 
chargee  of  the  lone,  I  ie  the  ionic  strength,  and  a  and  B  are 
constants  for  which  empirical  values  (frequently  unity)  are 
assigned  [25,31]. 


HEAT  CAPACITIES 

High  temperature  data,  especially  heat  capacities,  for 
aqueous  reactions  are  available  in  several  forms.  Various 
approaches  which  are  used  to  obtain  and  manipulate  these  data 
will  be  susuaarized,  so  that  all  data  can  be  in  the  same  form  for 
systematic  use.  The  equations  which  define  the  important 
relationships  of  equilibria  and  thermodynamio  functions  are: 

-HT  In  E(T)  -AG* 

40*  -iSf-TiS*- 

4  ®5».  '  4  S»«  41  Ml4  \  "  -j£4  % 4  “  T 

The  equilibrium  constant,  K,  at  temperature,  T,  can  thus  be 
calculated  with  the  free  energy  at  TA  A  G*.  ,  or  with  the  values 
for  enthalpy,  A  H*  ,  and  entropy,  A  ,  it  that  temperature. 

The  equilibrium  constant  can  also  be  calculated  with  the  more 
accessible  values  of  free  energy  at  298*K,  and  the  standard 
state  entropy  change  at  298*K,  if  the  value  for  the  change  in 

heat  capacity,  A  C  ,  is  known  for  the  temperatures  of  interest. 

P 

For  most  systems,  heat  capacity  is  variable  and  values 
are  least  accessible  at  high  temperatures.  Where  data  are 
available,  C  is  usually  defined  over  a  temperature  range 
(specific  topthe  compound  in  question)  by: 

AC  -  a  +  bT— 1  +  cT±2  (B-6) 

P 

As  these  data  are  frequently  not  available,  three  means  of  making 
estimates  are  used  [25,32].  These  are  the  zero  heat  capacity 
assumption,  the  constant  heat  capacity  method,  and  the 
correspondence  principle.  Over  the  wide  range  of  temperatures 
occurring  in  boilers  it  is  unlikely  that  AC*  will  be  zero. 
However,  if  two  chemical  reactions  are  combined  such  that  like 
charges  are  balanced  it  is  occassionally  possible  that  the 
changes  in  heat  capacity  will  cancel,  yielding  a  situation  where 
the  overall  change  in  heat  capacity  is  zero.  This  is  true  for 
the  ionization  of  dihydrogen  phosphate: 

HjPO,"  -  H+  +  HPo/~  K  (B-7) 

H+  +  0H~  -  H,0  1/^  (B-8) 


H,P04"  +  OH-  -  HP04  "+  H,0  K/Zw  (B-9) 


(B-3) 

(B-4) 

(B~5) 
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Snell  xero  bent  capacity  treatment  is  rarely  possible.  It  is  aueb 
■ore  likely  that  the  combination  of  two  chenical  reaetions  will 
lead  to  a  constant,  non-sero  heat  capacity.  In  this  constant 
heat  capacity  approach,  the  equations  are  coabined  so  that  the 
ionic  charges  are  balanced.  This  is  aore  likely  to  result  in 
constant  values  of  AC  if  the  species  are  siailar.  The  walidity 
of  this  approach  can  pbe  verified  by  plotting  known  values  and 
those  that  are  calculated  in  this  way.  There  is.  however,  no 
general  test  of  this  aethod.  It  is  thus  used  with  the  greatest 
confidence  to  extrapolate  values  to  higher  teaperatures  where 
soae  higher  teaperature  aeasureaents  have  been  aade. 

The  third  approaoh  to  estiaating  heat  capacity  values  is  the 
correspondence  principle  described  by  Criss  [32]  and  Cobble  [331. 
For  this  approach  a  heat  capacity  funotion  is  defined: 

C*  -  AT  *  S*  S*  dT  (B-10) 

where  S*.  is  defined  as: 

^  "  *T  +  bTS$98  (B"U) 

This  heat  capacity  term  can  then  be  used  to  approximate  standard 
free  energy  by: 

®r  -  53«-  “  *  V 

Evaluations  of  this  approach  [32]  have  deaonstrated  that  it  is 
satisfactory  up  to  250*C;  at  300*C  there  is  still  considerable 
uncertainty. 
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APPENDIX  C 


HIGH  TEMPERATURE  THERMOCHEMICAL  DATA 


To  model  boiler  water  it  ie  necessary  to  compile  and  assess 
tbe  available  high  temperature  thermochemical  data.  A  wide 
variety  of  sources  were  nsed  to  collect  these  data;  as  a  resnlt 
the  data  are  given  in  a  variety  of  forms.  For  example,  the 
United  States  Geological  Survey  (USGS)  [34]  compiled  data  on  the 
energy  for  the  reactions  of  formation  from  the  elements  of  many 
species  of  interest  here.  Values  for  AH® ,  S# ,  AG# ,  and  log  K.  at 
temperatures  ranging  from  298°K  up  to  1800°K  for  some  compounds 
are  given.  Additionally,  Helgeson  [35]  has  compiled  iot. 
dissociation  equilbria  for  aqueous  species,  solubility  products, 
heat  capacity  values,  and  enthalpies,  entropies,  and  average  heat 
capacities  for  the  formation  of  aqueous  ions;  data  are  given  up 
to  300° C  where  they  are  available. 

Data  are  also  included  in  reports  of  modeling  done  of  boiler 
chemistry  [10-12,31]  and  of  natural  waters  [22,38].  Most  of 
these  data  are  presented  as  values  of  the  formation  constant  or 
the  solubility  product  at  a  specific  temperature.  However,  some 
of  the  data  are  presented  in  the  more  general  form: 

log  K  -  a  +  bTiX  +  cT±2  (C-l) 

Three  further  sources  [25,37,38]  include  values  of  heat 
capacities  and  of  log  K  at  different  temperatures.  Finally,  some 
applicable  thermodynamic  data  are  available  in  the  form  of  phase 
diagrams  [3,5]  . 

The  data  found  in  these  sources  have  been  summarized  in 
tables  of  Appendix  C,  for  the  purposes  of  this  work,  by  species 
type.  As  water  is  the  medium  in  which  all  solution  reactions  in 
the  boiler  take  place,  data  for  the  ionization  of  water  are 
essential.  Many  sources  of  these  data  are  available  and  the 
reported  values  are  in  Table  C-l.  Data  for  other  species  follow 
in  several  tables;  to  make  these  tabulations  easier  to  use. 

Table  C-2  includes  common  general  information.  There  are  few 
data  pertinent  to  phosphoric  acid  and  phosphate  species,  although 
the  first  two  ionization  constants  have  been  carefully  measured 
[25,31].  Other  pertinent  data  are  in  Table  C-3 .  One  different 
source  of  data  is  available,  that  is,  phase  diagrams  for  sodium 
phosphate  [3,5].  A  study  of  the  available  data  leads  to  the 
conclusion  that  there  is  not  much  high  temperature  data  for  metal 
phosphate  species.  Because  of  the  possible  implication  of  iron 
phosphate  species  involvement  in  phosphate  hide-out  [5]  and 
because  of  the  known  interaction  of  phosphate  with  iron  in 
natural  water  systems  [39] ,  these  missing  data  may  be  important 
to  the  development  of  a  complete  model.  Consideration  of  various 
phosphate  forms  is  apparently  unnecessary  as  they  are  converted 
to  orthophosphates  at  high  temperatures  [40]  . 


Data  for  CO, ,  bicarbonate,  and  carbonate  reactions  are 
widely  available.  They  are  summarized  in  Table  C-4.  Table  C-5 
includes  data  for  the  reactions  of  SiO,,  silicic  acid,  silicate, 
and  aetal  silicates.  Data  for  hydroxide  species  and  metal 
hydrolysis  reactions  are  in  Table  C-6 .  Reactions  involving 
sulfuric  acid,  hydrogen  snlfate,  and  aetal  sulfates  have  been 
well  studied;  data  are  in  Table  C-7 .  Data  for  species  involving 
chloride  which  are  aainly  ion  pair  dissociation  reactions,  have 
been  presented  in  Table  C-8.  These  data  are  more  important  at 
high  temperatures  than  at  low  because  there  is  a  greater  tendency 
for  ion  pair  association  due  to  the  decreasing  solvating  capacity 
of  water.  The  energies  of  formation  of  metal  oxides  are 
available  [34,35]  for  FeO,  Fe,0, ,  Fe,04,  MgFe,04,  and  Cu,0.  The 
only  reported  solubility  data  for  metal  oxides  are  for  magnetite, 
Fe,04  [31] .  Enthalpies  of  formation,  with  entropy  and  average 
heat  capacities  are  available  for  the  formation  of  aqueous  metal 
ions  and  a  variety  of  anions  [35].  These  ions  are:  Ag+,  Cu+, 
Cu*+,  Pb*+,  Zn*+,  Fe»+.  Fe*+,  Ca*+,  Mg*+,  Ba*+,  K+,  Na+,  Al*  +  , 
HS-,  S*-,  F-,  S04*-,  CO,1-,  OH- ,  C1-,  and  H+.  These  data  are 
useful  in  deriving  other  values  for  aqueons  media.  Notably 
absent  in  this  source  are  data  for  F04*-.  Included  in  Table  C-9 
are  data  for  all  other  species  which  did  not  fit  the  above 
classifications. 
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TABLE  C-l 


DATA  ON  THE  IONIZATION  OF  WATER  (  -log  ZJ 
Reaction:  H.O  »  H+  +  OH- 


Teaperatnre  (#C)  pK^4 1  pK^44 

0  14.94 

25  13 .99 

50  13 .27 

60 

75  12.71 

100  12.26 

125  11.91  11.942 

150  11.64  11.675 

175  11 .44  11 .477 

200  11.30  11.352 

225  11.22  11.261 

250  11.20  11.230 

275  11.22  11.245 

300  11.30  11.300 

325 

350  11.521* 

352 

374.15  11.674* 

*  Extrapolated  beyond  the  range 


PV* 

PV* 

’V* 

>V‘ 

11. 

14.938 

14.00 

13.995 

13.27 

13.275 

13.03 

12.712 

12.26 

12.265 

11.907 

11.905 

11.912 

11.634 

11.64 

11.671 

11.638 

11.420 

11 .461 

11.432 

11.254 

11.27 

11.302 

11 .289 

11 .131 

11.196 

11.208 

11.049 

11.13 

11.138 

11 .191 

11.012 

11.132 

11 .251 

11.034 

11.39 

11.193 

11 .406 

11.71 

11.42 

12.30 

11.15 

15.74 

11.997* 

of  the  data. 


log  ([Hf][OB-])  -  -  151.713/T  -  111.491  -  0.03685  T  + 

44.077  log  T  +  (2  A  •  I1/2/  (  1  +  1.5*I1/2)) 
-(  0.6356  -  0.001078)  I  (from  reference  12) 


TABLE  C-2 


SUMMARY  INFORMATION  FOR  THE  THEKMODTNAMIC  DATA  TABLES 


All  species  ire  assumed  to  be  dissolved  aqueous  species  unless 
otherwise  noted. 

*  indicates  that  data  at  other  temperature  values  are  available 
in  the  reference.  Usually  the  value  at  the  highest  temperature 
is  presented. 

a(T)  -  AT  +  Tt  ln(  T,  /  T4) 

•  in  a  chemical  formula  signifies  a  dissolved  ion  pair. 

Dimensions  are:  AH,  cal/mole;  A6,  cal/mole;  AS,  cal/aole-deg; 

C  ,  cal/mole-deg. 

P 

Formation  data  available  are:  AH,  S° ,  AG,  C*,  and  log  K.  for 
reactions  of  formation  from  the  elements  inptheir  satanaard 
states. 

Temperatures  used  are  all  *Cj  thus  log  E(200)  signifies  the 
formation  constant  at  200SC. 


Subscript  s  ”  solid 


TABLE  C-3 


THERMO CHEMICAL  DATA  FOR  PHOSPHORIC  ACID  AND  PHOSPHATES 


Reference 

Reaction 

Data 

23 

H+  +  P04‘-  -  HP04*- 

A  H,  ■  -3 .50  A  Sa  “  43 

23 

fff  +  hpo4»-  -  h,po4- 

A  Rt,  -  -0.80  A  Saa  -  30.3 

23 

Hf  +  H,P04-  -  H,P04 

A  Hla  «  1.88  ASaI  -  16.0 

25,31 

H,P04  -  Hf  +  H4P04- 

pK  (300)  -  4.400* 

23,31 

HjP04-  -  Hf  +  HP04»- 

pi  (300)  -  8.946* 

25 

QH-  +  HaP04-  * 

AGa>,  -  -9271  ASa,,  -  -11.0 

HP04»-  +  HaO 

AC  (0-300*0  -  -0.1 

25 

H,P04  +  OH-  = 

AGaaa  -  -16163  ASat,  *=  -3 

HaP04-  +  Ha0 

AC  (0-300*0  -  11 

34 

h,po4 

formation  data  to  700*C 

31 

NaaP04 

solubilty  1500  ppn  (350*0 

2,4 

sodium  phosphate 

phase  diagrams 

34 

Ca»(P04)» 

formation  data  to  1300*C 

34 

Ca,(P04),OH 

formation  data  to  1200*C 

26 

Fe*+  +  HP04»-  -  FeHP04 

log  8(30)  -  7.3 

26 

MgJ+  +  HP04»-  -  MgHP04 

log  0(25)  -  2.91 

26 

Ca*+  +  HP04*-  -  CaHP04 

log  0(25)  -  2.70 

26 

Ca*+  +  HaP04»-  -  CaHaP04+ 

log  0(25)  -  1.08 

26 

Cn*+  +  HP04*-  -  CnHP04 

log  0(25)  -  3.2 

26 

Na+  +  HP04*-  «  N«HP04- 

log  0(25)  -  0.60 

26 

K+  +  HP04»-  -  OP04- 

log  0(25)  -  0.49 

TABLE  C-4 


THERMOCHEMICAL  DATA  FOR  00, .  HCO,-,  mod  CO,*-  SPECIES 


Reference  Reaction 

Data 

25 

w*(a) 

«>*;:;  -  «*(.,> 

C*  (25-300*0)  -  45 

34 

formation  data  to  1500°C 

25 

AO  -  1992  +  23.0  AT 
+  36.0  a(T) 

25 

00,  +  OH-  -  HCO,- 

AO  -  -10406  +  3 .7  AT 

-  19.5  a(T) 

25 

H00,-  +  OR-  * 

AO  -  -5010  +  16.1AT 

00,*-  +  H,0 

-  21.0  a(T) 

22 

00,*-  +  »■  »  HOO,- 

log  X  -  -  6.4980  +  0.02379 

+  2902. 3  9/T 
(valid  0-50*0 

11 

HOO,-  -  BH-  +  00,*- 

log  1(250)  -  -11.34 
log  E(300)  -  -11.98 

35 

HOO,-  -  *•  +  00,*- 

log  E(250)  -  -11.43 
log  X(300)  -  -13.38 

22 

HCO,-  +  H+  -  H,CO, 

log  I  -  -14.8435  +  0.03279 

+  3403 .71/T 
(valid  0-50*0 

11 

H,CO,  -  B*  +  HOO,- 

pE(250)  -7.88  pC(300)  -8.58 

35 

H,CO,  -  H+  +  HOO,- 

pX(250)  -7.63  pX(300)  -8.86 

10 

H,CO,  -  JH-  +  HOO,- 

log  1(300)  -  -8.55 

12 

H,CO,  -  Bf  +  HCO,- 

log  E  -  -2382  .2/T  +  8.153 

-  0.02194  T 

35 

a)»*"(aq) 

‘‘C’e.lcit.)  " 

Ca*+  +  00,*- 

formation  data  up  to  300*C 

34 

formation  data  op  to  900*0 

22. 

,36 

log  E  -  13.543  -0.0401  T 
(va??8°6^50*C) 

36 

(s.argonite) 

Ca*+  +  00,*- 

log  E  -  -10.21  +  0.0217  T 
-5.17  *  10-*  T* 

35 

CaCO,(  -  Ca*+  +  00,*- 
CaOO,J#j  -  Ca*+  +  00,*- 

log  E(300)  -  -14.10 

11 

log  E(250)  -  -12.72 
log  X(300)  -  -14.10 

22 

Ca»+  ♦  CO,*-  -  CaCO,* 

log  E  -  -27.393  +  4114/T 

+  0.05617  T 
(valid  10-50*0) 

35 

CaOO,  -  Ca*+  +  00,*- 

log  K(200)  -  -5.2 

35 

CaCO,  +  2  »+■  - 

log  E(300)  -  8.73 

Ca*+  +  H,0  +  00, 

22 

Ca»+  +  HOO,-  -  CaHCO,+ 

log  E  -  -  2.95  +  0.0133  T 

(valid  0-50*0 

25 

Ca»+  +  HCO,-  -  CaHCO,+ 

pX(250)  --2.94 
pE(300)  -  -3.43 

(« 

if 

CaHOO,  i*  delated  thia  value 

will  change) 

TABLE  C-4  (CONT'D) 


THERMOCHEMICAL  DATA  FOR  CO, ,  HOO,-,  tad  CO,*- 


Reference 

Reaction 

Data 

34 

CaMg(CO,),(  , 
"•“.(a) 

formation  data  to  800*C 

34 

formation  data  to  700*C 

22 

Mg»+  +  HCO,- 

-  MgHCO,+ 

log  K  -  2.319  -  0.011056 
+  2.298  x  10-*  T* 
(valid  10-90*0 

25 

Mg*+  +  HCO,- 

-  MgHCO,+ 

pK(250)  -  -2.56 
pC(300)  -  -2.97 

22 

Mg*+  +  CO,»- 

-  MgCO, 0 

log  K  -  0.991  +  0.0066  T 
(valid  10-90*0 

25 

M*+  +  HOO,-  - 

MHCO,+ 

log  E  at  variona  T  for  Ni 
Mn,  Cu,  Zn,  Fe 

TABLE  C-5 


THERMO CHEMICAL  DATA  FOR  SILICA, 
SILICIC  ACID,  AND  SILICATE  SPECIES 


l«f*ren<i 

Reaction 

Data 

25 

SiO,.  .  +  2  H,0  - 
1,1  H4Si04 

AGC298)  -  5494  AS(298)  -5.26 

AC  (0-300)  -  -11.9 

35 

SiO„  .  +  2  H,0  « 

' '  H4S104 

logPE(300)  -  -1.94* 

25 

SiO>(  ,  +  2  HaO  - 

H*Si0*uq) 

AG  -  5494  -  5.26  AT 

-11.9  a  (T) 

pE(300)  -2.00  • 

25 

SiOa,  .  +  2  H,0  - 

H,Si04-(  , 

H4Si04  +  OH-  - 

pE(300)  -  11.22* 

25 

AG(298)  -  -5687  AS(298)-5.14 

H,SiO,-  +  HaO 

AC  (0-300)  -  12 

25 

H4Si04  +  OH-  - 

logPr  -  2346 .69/T  + 

HjSiO,-  +  HjO 

2.47979  In  T  - 
18.4014 

25 

H4Si04  «  »■  +  HjSi04- 

pE(300)  -  9.22* 

22 

H4Si04  -  W-  +  H,Si04- 

log  E  -  6.3680  -  0.0163  T 

-  3405.8999/T 

22 

H4Si04*  -2  Hf  +  H,Si04»- 

log  E  -  39.4780  -  0.0659  T 

-  12355 .0977/T 

35 

H4Si04 

Conation  data  to  300*C 

25 

M1Si°4(,)  -MO(i)  + 

MSi°>(  , 

M  -  Mg 

pE(300)  «  2.182* 

K  -  Ca 

pE(300)  -  4.265* 

M  -  Fa 

pE(300)  -  0.426* 

25 

MSi0*(.)  "  K0(.)  + 

X  -  Xg 

pE(300)  -  3.242* 

K  -  Ca 

pE(300)  -  8.115* 

X  -  Fa 

pE(300)  -  2.068* 

25 

XSiOt(#)  +  HaO  -  M*+  + 

Si0*(.,qn.rt*)  +  2  <*' 

* 

M  -  Mg 

pE(300)  -  16.5* 

M  -  Ca 

pE(300)  -  14.2* 

M  -  Fa 

pE(300)  -  19.1* 

25 

M*Si0«(a)  +  H*°  “  2M*+ 

+  Sl0»(.,q»artz)  +  4  0H' 

- 

X  -  Mg 

pE(300)  -  31.9* 

X  -  Ca 

pE(300)  -  24.6* 

X  -  Fa 

pE(300)  -  36.7* 

35  +  4  B+  -  log  1(300)  -  12.61* 


2  Mg»+  +  H4Si04 


28 


TABLE  C-5  (CONT'D) 


i 

i 

- 

i 


* 


% 


THERMOCHEMICAL  DATA  FOR  SILICA. 
SILICIC  ACID,  AND  SILICATE  SPECIES 


•  1 


25  MSiO,(#)  +  H40  +  Hf  - 


M*+  +  H.Si04- 

M  -  Mg 

pl(300)  -  5.112* 

M  -  Ca 

pE(300)  -  2.854* 

M  -  Fe 

pl(300)  -  7.786* 

i 

i 

25 

M»+  +  N,Si04(#)  - 

N»+  +  MNSi04 »  . 

M  «  Fe,  N  -  Ci  ' 

pi (3 00)  -  -6.949* 

M  -  Ce,  N  -  Fe 

pl(300)  -  5.102* 

M  -  Mg.  N  -  Ce 

pE(300)  -  -5.211* 

- 

- 

M  -  Ci,  N  -  Mg 

pl(300)  -  2.026* 

P 

i 

34 

CiMg(SiO,)s  (#) 

CiSiO  } 

Cisioj*;, 

fomtioa  diti  to  1300*C 

34 

fomtioa  diti  to  1100°C 

34 

fomtioa  diti  to  1500°C 

34 

Ci,Mg(Si04),  (  j 

FeaSi04  (  , 

MgSiO 

C-.Si.oJJ^  +  7  HaO  = 

fomtioa  diti  to  1300*C 

34 

fomtioa  diti  to  1500°C 

i 

* 

34 

fomtioa  diti  to  1300*C 

34 

fomtioa  dati  to  1500*C 

12 

log 

R{275)  -  -43.16 

3Ci»+  +  6  GB-  +  2  H4Si04* 

12 

CiSiO, (#)  +  3  HaO  -  Ci»+ 

log 

K(275>  -  -15.70 

r 

( 

-  -  — 

+  2  OH-  +  H4Si04* 

12 

CiaSi04<sj  +  4  HaO  - 
2  Ci*+  +  4  OH-  +  H4Si04* 

log 

1(275)  -  -25.33 

12 

Ci,SiOf^#j  +  5  HaO  - 
3  Ci»+  +  6  OH-  +  H4Si04* 

log 

K(275)  -  -29.6 

► 

_ « 

12 

***si°.(,)  +  3  H,0  *  Mg*+ 

log 

1(275)  -  -18 

+  2  OH-  +  H4Si04* 

12 

3*g»Si»07 .2HaO,  .  +  5  HjO 
-  3  Mg*+*+  6  oh- 

log 

1(275)  -  -52.14 

; 

+  H4Si04° 

i 

• 

12 

M8*Si04(  .  +  *  H,0  » 

2  Mg»+  +  4  OH- 

log 

1(275)  -  -32.77 

+  2  H4Si04# 

12 

CiMgSia04(#)  +  6  HaO  - 

log 

1(275)  -  -35.47 

Ci»+  +  Mg»+  +  4  0H- 
+  2  H4Si04# 

i 

f 

12 

CijMgSi jO,^s^  +  7  HaO  - 

log 

K275)  -  -45.46 

i 

2  Ci*+  +  Mg*+  +  6  OH- 
+  2  H4Si04* 


12 

CiMgSi04(#)  +  4  H,0  - 

log  1(275)  -  -30.68 

Ci*+  +  Mg»+  +  4  OH- 

+  H4Si04. 
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*  J 


TABLE  C-6 


THERMOCHEMICAL  DATA  FOR  HYDROXIDES 
AND  METAL  HYDROLYSIS  REACTIONS 


Reference 

Reection 

Data 

35 

0H-,  . 

(aq) 

formation  data  to  300°C 

31,12 

Ca(OH)l(s)  -Ca*+  +  2 

0H- 

log  Kip  -  -25.7085 

+  12.3722  log  T 

-  530.49/T 

-  0.032331  T 

+  (6A  I1/2/l  +  Ak  I1/2) 

-  BI  -  Cl*  (see  26  for 
parameters,  valid  to 
oritical  T  of  water) 

11 

Ce (OH) a %  =Ca»+ 

'*'+  2  oh- 

pK(250)  -  8.37 
pE(300)  =  9.38 

25 

.  -Ca»+  +  2 
Ca(OH)J*{  -Ca»+  +  2 
Ca(OH)aJ  J  -Ca»+  +  2 

OH- 

pE(300)  -  10.48* 

10 

OH- 

log  X(300)  -  -9.39 

12 

OH- 

log  K(275)  -  -8.88 

25 

Ca*+  +  2  H#0  » 

Ca  (OH) . .  »  +  2  H+ 

pK(300)  -  12.12* 

34 

Ce(OH) 

Mg(OH>a  ; 

Mg(OH)aJ‘;  -Mg*+  +  2 
Mg(OH),J*j  -Mg*+  +  2 

formation  data  to  400° C 

34 

formation  data  to  600°C 

25 

OH- 

pK(300)  -  14.54* 

11 

OH- 

pl(250)  -  13.97 
pE(300)  -  14.98 

31 

Mg(OH)a ,  .  -Mg*+  +  2 
Mg(OH),J‘:  -Mg*+  +  2 
Mg(OH)a}*J  -Mg*+  +  2 
Mg(OH),J*j  -Mg»+  +  2 

0H- 

pK(300)  -  14.98* 

10 

OH- 

log  1(300)  -  -14.25 

12 

OH- 

log  r  -log  [Ca*+] [OH-] *  -5.6 

12 

OH- 

log  E(275)  -  -14.48 

35 

Mg (OH) +  -  Mg*+  +  OH- 

pK(200)  -  4.1* 

10 

Mg (OH) +  -  Mg*+  +  0H- 

pK(250)  -  5.0  pX(300)  -  5.7 

22 

Mg*+  +  OH-  «  Mg(OH)+ 

log  X  -  0.684  +  0.0051  T 
(valid  10-90*0 

25 

Mg*+  +  2  H,0  « 

Mg (OH) .  (s)  +2  Hf 

pK(300)  -  8.06* 

12 

NH40H  -  nh4+  +  oe- 

log  K  -  -2.8  x  10-*  T* 

+  0.018  T  -  7.63 

+  A  I1/2/  (1  +1.5  I1/2) 

31 

NH40H  -  NH4+  +  OE- 

pl_(300)  -  6.466* 

10 

NH40H  -  NH4+  +  OH- 

log  E(300)  -  -6.21 

30 


TABLE  C-6  (CONT'D) 


THERMO CHEMICAL  DATA  FOR  HYDROXIDES 
AND  METAL  HYDROLYSIS  REACTIONS 

Reference  Reaction  Data 


31 

NaOH  -  Na+  +  OH- 

K(218)  -  0.55 

10 

NaOH  •  Na+  +  OB~ 

log  K(300)  -  -1.10 

31 

LiOH  -  Li+  +  OH- 

1(350)  -  0.007 

..  - 

35 

FeOH+  -  Fe*+  +  QH- 

pK(200)  -  5.7* 

• 

25 

Fe*+  +  H.O  -  Fe(OH)+ 

AG(298)  «  12700 

+  Hf 

AS(298)  -1.7  AC  (0-300)  —7 

P 

42 

2  Cn*+  +  2  HaO  - 

log  X  -  2.497  -  3833/T 

Cn,(OH)»*+  +  2  Hf 

25 

M*+  +  HaO  -  M(OH)+  +  Hf 

pE(300) 

• 

M  -  Mn 

5.77* 

M  -  Fe 

4.92* 

M  -  Co 

5.10* 

M  -  Ni 

5.22* 

M  »  Cu 

4.22* 

25 

M»+  +  2  HaO  -  M(OH)a*  +  2 

Hf  pK(300) 

9 

M  -  Mn 

11.9* 

H  -  Fe 

10.9* 

M  -  Co 

10.2* 

M  -  Ni 

11.6* 

M  -  Cu 

9.2* 

L  _ 

25 

M*+  +  3  H.0  «  M(0H) a-  +  3 

Hf  pC(300) 

9 

M  -  Mn 

21.0* 

M  -  Fe 

18.2* 

M  -  Co 

19.3* 

M  «  Ni 

18.1* 

M  -  Cu 

17 .4* 

,, 

9 


31 


t 


TABLE  C-7 


THERMO  CHEMICAL  DATA  FOR  SULFURIC  ACID,  AND  SULFATE  SPEaES 


Reference 

Reaction 

Data 

35 

S04»- 

foraation  data  to  300*C 

25 

hso4-  «  B+  +  S04*- 

pK(300)  *  6 .324* 

10 

HS04-  -  fff  +  S04»- 

log  1(300)  -  -7.06 

11 

HS04-  -  fff  +  S04»- 

pK(250)-  4.98  pK(300)-  5.71 

35 

hso4-  -  H+  +  S04»- 

log  R(300)  -  -7.06* 

31,12 

hso4-  -  H+  +  S04x- 

log  E  -  -  33 .0024  log  T  + 
91.471  -  3520 .3/T 

22 

H+  +  S04»-  -  HS04- 

log  K  -  -5.3505  +  557.2461/T 
+  0.0183412  T 
(valid  25-225*0 

25 

S04»-  +  HtO  - 

HS04-  +  OH- 

pK(300)  •  4.977* 

10 

H,S04  «  H+-  +  HS04- 

log  1(300)  -  0.70 

35 

H#S04  ■  fff  +  HS04- 

log  K(300)  -  0.50* 

25 

hso4-  +  hpo4»-  - 

h,po4-  +  so4»- 

pK(300)  -  2.622* 

31,12 

CaS04(#)  -  Ca*+  +  S04»- 

log  K  -  -  133.207  +  3569.6/T 
+  53.5472  log  T 
+  0.0529025  T 

25 

CaS04f  ,  -  Ca*+  +  S04»- 
CaS04[  (  -  Ca*+  +  S04*- 
CaS04  «  Ca*+  +  S04*- 

pl(300)  -  5.917* 

10 

log  K(300)  -  -  9.05 

10 

log  K(300)  -  -  5.36 

35 

CaS04  -  Ca*+  +  S04*- 

log  K(200>  -  -  3.6* 

12 

CaS04f  .  -  Ca»+  +  S04»- 
CaS04J*j  -  Ca*+  +  S04*- 
B.S04  ,  -  Ba*+  +  S04*- 
BaS04J*j  +  Ca*+  - 

17  CaS04f  .  +  Ba*+ 
MgS04*  -  Mg*+  +  S04»- 

log  K(275)  «  -  9.03 

11 

pf(250)  -8.03  pK(300)  -9.25 

25 

pE(300)  »  10.155* 

25 

pK(300)  -  4.238* 

31 

log  K  -  -158.540  +  4180.6/T 
+  62.160  log  T 
-  0.046298  T 

10 

MgS04  -  Mg*+  +  S04»- 

log  K(300)  -  -7.92 

11 

MgS04 ( s )  -  Mg*+  +  S04»- 

pK(250)  -5.56  pK(300)  -7.33 

35 

*gS04  -  Mg*+  +  S04»- 

log  E(200)  -  -4.8* 

31 

Na+  +  S04>-  -  NaS04- 

pK(300)  -  2.84* 

10 

NaS04-  -  Na+  +  S04*- 

log  K(300)  -  -  3.64 

10 

NaaS04  «  2  N*+  +  S04*- 

log  K(300)  -  2.40 

11 

NaaS04  »  2  Na+  +  S04*- 

pK(250)  -5.01  pE(300)  -5.34 

10 

NaHS04  -  Na+  +  HS04- 

log  1(300)  -  -0.30 

22 

z+  +  so4»-  -  rso4- 

log  r  -  3.1060  -  613 .5999/1 
(valid  10-50*0 

31 

1+  +  so4»-  «  kso4- 

pK(250)  -  2.45 

35 

kso4-  «  k+  +  so4»- 

log  1(200)  -  -1.94* 

35 

ihso4  -  r+  +  hso4- 

log  1(250)  -0.8 
log  E(300)  -  -  0.3 

10 

nh4so4-  -  nh4+  +  so4»- 

log  E(300)  -  -0.30 

35 

MnS04  -  Mn*+  +  S04>- 

log  K(200)  -  -4.3* 

35 

ZaS04  -  Zn*+  +  S04»_ 

log  E(200)  -  -4.6* 

TABLE  C-8 


THERMO CHEMICAL  DATA  OF  CHLORIDE  SPECIES 


Reference 

Reaction 

Data 

• 

4 

35 

Sci*’1 

FeCl,  J 

FeCl,  J 

HC1(.) 

formation  data  up  to 

300*C 

34 

formation  data  np  to 

700*  C 

34 

formation  data  up  to 

1000* C 

34 

formation  data  np  to 

400*C 

_ _ _ 

_ „ 

34 

formation  data  np  to 

1500*C 

4 

34 

NH4C1(  J 

HC1  «  H+  +  Cl- 

formation  data  np  to 

400°C 

10 

log  K(300)  -  -1.23 

35 

HC1  *  H+  +  Cl- 

log  K(300)  -  -1.24* 

31 

HC1  -  H+-  +  Cl- 

log  K(360)  «=  -2.98* 

■ 

10 

NaCl  -  Na+  +  Cl- 

log  K(300)  -  -1.10 

• 

< 

4 

35 

NaCl  -  Na+  +  Cl- 

log  K(300)  -  -0.8* 

31 

NaCl  -  Na+  +  Cl- 

log  K(360)  -  -1.88* 

11 

NaCl .  .  -  Na+  +  Cl- 
CaCljJ'j  -  Ca*+  +  2C1- 

log  K(250)  -.9  log  K(300)*.3 

11 

log  K(250)  -  2.4 
log  r(300)  -  1.6 

11 

MgC1*(«)  “  M**+  +  2C1~ 

log  1(250)  -  1.4 
log  1(300)  -  0.6 

:• 

35 

KC1  «  R+  +  Cl- 

log  E(300)  -  -0.6* 

35 

FeCl*+  «  Fe»+  +  Cl- 

log  0(150)  -  -3.98* 

35 

FeCla+  «  Fe*+  +  2C1- 

log  0(150)  =  -4.72* 

35 

FeClj  -  Fe*+  +  3C1- 

log  0(150)  -  -4.30* 

35 

FeCl4-  «  Fe*+  +  4C1- 

log  0(150)  -  -3.23* 

-•  ~~  " 

* 

•  .  m 
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3 


TABLE  C-9 


Reference 

34 

2S 

25 

22 


35 

34 

34 

25 


10 

11 

22 


25 

25 

25 

25 

25 

25 


25 

25 


THERMO  CHEMICAL  DATA  FOR  OTHER  SPECIES 


Reaction 


Data 


°*<aq) 

°*<g> 

H,S  - 


"  °*(aq> 
Hf  +  HS- 


HS-  -  H+  +  S* 


H,S  - 
NHi(  , 

m>\\) 


H+  +  HS- 
“  (aq) 


"“•(aq)  "  "“‘(g) 
NH4+  «  NH,  +  H4- 


NH4+  »  NH, •  +  H+ 


"“•(aq) 

NH,  +  H,0  -  NH4+  +  OH- 
NH,  +  Hf  «  NH4+ 

NH,(b)  +  H,0  -  NH4+  +  0H- 

“*(aq) 

“*(g)  "  “*(aq) 

"•(aq) 

"•(g) N*(aq) 


formation  data  to  1500°C 

C  (25-300)  -  45 
P 

AG-  -  3944  +  22.9  AT 
+39  a(T) 

log  K  -  11.17  -  0.02386T 
-3279. 0/T 
( valid  25-300*0 
log  K(300)  -  -7.72* 
log  K(300)  -  -8.06* 
formation  data  to  1500°C 
AG-  -  -2410  +  19.4  A  T 
+5.9  a(T) 
log  D(300)  -0.7 

log  K(250)  -  -5.14 
log  1(300)  m  -4.57 
log  E  -  0.6322  -  2835 .7598/T 
-  0.001225  T 

C (25-300)  »  14 
P 

pK(300)  -  6.694* 

pK(300)  -  -4.607* 

pK(300)  »  7.582* 

C  -  36 
P 

AGT  -  4238  +  17 .5  AT 
+  29  a(T) 

C  -  46 
P 

AG-  -  4340  +  23 .0  AT 
+39  a(T) 
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